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ABSTRACT 

We present models for the slow neutron-capture process [s process) in asymp- 
totic giant branch (AGB) stars of metallicity [Fe/H]=— 2.3 and masses 0.9 Mq 
to 6 Mq. We encountered different regimes of neutron-capture nucleosynthesis 
increasing in importance as the stellar mass decreases: the ^^Ne(a,n)^^Mg reac- 
tion activated during the thermal pulses, the ^^C(a,n)^^0 reaction activated in 
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radiative conditions during the interpulse periods, and the ^^C(a,n)^^0 reaction 
activated during the thermal pulses, also as a result of mild proton ingestion 
episodes. The models where the ^^C burns radiatively (masses ~ 2 Mq) produce 
an overall good match to carbon-enhanced metal-poor (CEMP) stars showing s- 
process enhancements (CEMP-s), except they produce too much Na and F. On 
the other hand, none of our models can provide a match to the composition of 
CEMP stars also showing rapid-process enhancements (CEMP-s/r). The models 
fail to reproduce the observed Eu abundances, and they also fail to reproduce 
the correlation between the Eu and Ba abundances. They also cannot match 
the ratio of heavy to light s-process elements observed in many CEMP-s/r stars, 
which can be more than ten times higher than in the solar system. To explain 
the composition of CEMP-s /r stars we need to invoke the existence of a different 
"s/r" neutron-capture process either with features in-between the s and the r 
processes, or generated by superpositions of different neutron-capture processes 
in the same astrophysical site or in sites linked to each other - for example, in 
multiple stellar systems. 

Subject headings: Nucleosynthesis, Abundances, Stars: Abundances, Stars: AGB 
and Post-AGB 



Introduction 



Following on from the discovery of a large number of halo stars with peculiar chemical 
compositions the last decade has seen much interest sparked in detailed spectroscopic ob- 



servations and nuc leosynthetic models of stars of low metallicity ( iBeers fc Christliebl 12005 



Sneden et al.l 120081 ). These peculiar compositions provide important constraints on the evo- 
lution and nucleosynthesis of the low-metallicity stars that existed in the early Galaxy. By 
interpreting the observations in the light of stellar nucleosynthesis models we can derive clues 
as to the sites and nature of the first nucleosynthetic processes in the Galaxy. This may ul- 
timately provide important insights into the form of the initial mass function of the first 
stars. The composition of metal-poor stars may also provide constraints on cold dark mat- 
ter cosmological models of the formation of galaxies, which predict that galac tic halos may 
be co mposed of different components of possibly different composition (e.g. 
2007h . 



CaroUo et al. 



Roughly 10-20% (ICohen et al.ll2005l : lLucatello et al.ll2006l : ICarollo et al.ll201l[ ) of old halo 
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stars are enriched in carbon ([C/Fejl]> 1) and are referred to as carbon-enhanced metal-poor 
(CEMP) stars. Most CEMP stars are also enriched in nitrogen, with [N/Fe]>l. In most 
cases [N/Fe]<[C/Fe], but also [N/Fe]~[C/Fe] and [N/Fe]>[C/Fe] are found in roughly 10% 
and 5%, respectively, of the total number of CEMP stars. Roughly 2/3 of CEMP stars 
also show enrichm ents in the elements heavier than iro n, which are produced by neutron- 
capture processes (lAoki et al.l 120071 : ISneden et al.ll2008l ). These enhancements range from 
pure s enhancements, of elements made primarily via s/ow-neutron captures (the s process) 
in the case of "CEMP-s" stars ([Ba/Fe]>l), to s and r enhancements, of elements made via 
the s process and the r process, i.e., rapid neutron captures (the r pr ocess) in the case of 
"CEMP-s/r" stars with [Eu/Fe]>l and [Ba/Eu]>0 Jjonsell et al.ll2006h . There are also halo 
stars with pure r enhancements, so-called rl (usually classifie d by setting 0.5<[E u/Fe]<l) 
and rll ([Eu/Fe]>l), one of which, the famous CS 22892-052 JSneden et al.ll2003h . belongs 
to the CEMP star family. CEMP stars without enhancements in the elements heavier that 
Fe are commonly referred to as "CEMP-no" stars. 

To deal with this wide variety of compositions different scenarios for the origin of halo 
stars with peculiar compositions have been proposed. Halo stars showing pure r enhance- 
ments are believed to have been born from material heavily polluted by a primordial core- 
collapse supernova (SN of Type II, SNII) where the r-process elements are supposed to have 
formed. Though we note that there are still no SNII models able to reproduce the conditions 
needed for the r process to occur. One of the most favored sites to date, the neutrino wind 
from the nascent neutron star, has rec ently been shown to have major difficulties in producing 
r-process nuclei ([Roberts et al.ll2010l ). Qualitatively, the C and s enhancements in CEMP-s 
and CEMP-s /r stars can be explained by the presence of an originally more massive binary 
companion to the observed star, which produced C and s-process elements during its asymp- 
totic giant branch (AGB) phase and transferre d them to its compani on via mass transfer, by 
either Roche-Lobe overflow or wind accretion (ILucatello et al.ll2005l ). About 50% of CEMP 
stars showing an s-process signature also show an r-process signature and are thus classified 
as CEMP-s/r stars. It is puzzling how CEMP-s/r stars could have formed in such a large 
number given that we believe the s and the r processes to be independent events occurring 
in different astrophysical site s, AGB stars and SNII, respectively. M any scenarios have been 
proposed (as summarised by lJonsell et al.ll2006l : iLugaro et al.ll2009l ) and most of them seem 
unlikely. Predictions from a scenario where the r-process enhanceme nts are suppos e d to be 
primordial to the binary system have been investigated in detail by iBisterzo et al.l (j2011bl ) 
and will be discussed in § 14. 1[ 



^The standard spectroscopic notation [X/Y]=logio(X/Y)— logio(X/Y)0 is used throughout this paper. 
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The N enhancements in CEMP stars are also mysterious. The highest enhancements 
[N/C]>0, may be due to the effect of proton captures at the base of the con yective envelope 
in AGB stars with ma sses higher than ~ 3 Mq (hot bottom burning, HBB lLattanzidll991 



Boothroyd et al.lll993l ). which very efficiently convert C into N. The milder enhancements 



can be qualitatively explained by H and He burning occurrin g together when protons are 



ingested in the He-b urning convective regions of AGB stars (jCampbell fc Lattanzid 12008 



Iwamoto et al.l |200J) , or by the operation of some form of deep mixing of the envelope 



material of the AGB star companion down to regions where H burning occurs. No physical 
mechanism for this de ep mixing has b een shown to be capable of producing the required 
nitrogen enhancement. IStancliffd (120101 ) has shown that it cannot be caused by thermohaline 
mixing. Another issue that has received much attention is the question of whether material 
transferred from the AGB companion onto the star now observed as a CEMP star was diluted 
on the secondary star via thermohaline mixing (or any non-convective process), or kept at 
the surfa ce until the star reached the giant branch and de veloped an extended convective 
envelope Jstancliffe et allboO?! : Istancliffe k Glebbeeklboosh . 

In this context, it is obvious that detailed models of the s process in AGB stars of low 
metallicity are crucial to improving our under st anding of CEMP stars . So me models have 



been p resented so far by IVan Eck et al.l (120031 ) , ICristallo et al.l (l2009b[). a n dlBisterzo et al. 



(I2OIOI ) and have been compared to stellar observations by iBisterzo et al.l ( l2011bl Jal). With 
this paper we greatly extend the range of published models of the s process in AGB stars by 
presenting the evolution and nucleosynthesis in a large set of stellar models of masses ranging 
from 0.9 Mq to 6 Mq at the typical metallicity of CEMP stars [Fe/H]=— 2.3. The s-process 
predictions for low- metallicity AGB stars of masses lower than 1.3 Mq and higher than 2 Mq 
are provided here for the first time in the literature. To test different physical and numeri- 
cal approaches we computed stellar evolutionary sequences using two different evolutionary 
codes. We then input these different evolutionary sequences into the same post-processing 
code, which employs a nuclear network from H to Bi and calculates nucleosynthesis by cou- 
pling mixing and burning in convective regions. We then compare our model predictions to 
the abundances of CEMP stars as a popu lation using the recent compilation of CEMP star 
compositions from iMasseron et al.l (l2010l ) . 



1.1. The s process in AGB stars of low metallicity 

From centre to surface an AGB star consists of: a C-0 degenerate core, a He-burning 
shell, a He-rich intershell, a H-burning shell, and an extended convective envelope, which is 
constantly eroded by strong stellar winds. When the whole envelope is lost the hot stellar 
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core is left to evolve into the post-AGB phase, then it becomes the central star of a planetary 
nebula and event ually a cooling white dwarf. For a detailed review on AGB evolution we 
refer the reader to iHerwigl (120051 ) . A peculiarity of the AGB evolution is that the H- and the 
He-burning shells are activated alternately. Most of the time the H-burning shell is active 
and freshly made He increases the mass of the intershell for periods of time of around 10^- 
10^ yr, depending on the stellar mass, with shorter periods associated with more massive 
AGB stars. The increasing mass of the intershell leads to an increase in the temperature 
at its base, until He burning is suddenly triggered resulting in a thermal pulse (TP). The 
energy released in this first phase of He-shell burning cannot be carried out radiatively and 
a convective region develops in the intershell, which lasts a few hundreds years (the exact 
duration depends on mass) and normally extends up to just below the H-burning shell. In 
low-mass and low-met allicity AGB stars this convective regi on can extend into the H- ric h 
layer and peculiar nucleosynthesis is expected to occur the re ( ICampbell fc Lattanzidl2008l ). 
including s- process nucleosynthesis ( ICristallo et al.ll2009al ). Mild proton-ingestion events of 
this type are found in our low-mass models. They have a strong impact on the s-process 
distributions and will be discussed in detail in § [31 

A TP causes the star to expand, which in turn causes the H-burning shell to be quenched. 
At this stage, the convective envelope can penetrate into the underlying layers and carry to 
the stellar surface material processed by H and He burning (a process referred to as the 
third dredge-up, or TDU for short), while He burning carries on radiatively. Eventually, He 
burning is also quenched and the whole cycle of H burning, TP, and TDU starts again. This 
cycle is repeated a few to hundreds of times depending on the time it takes for the mass loss 
to erode the envelope, which in turn depends on the initial stellar mass, and possibly the 
metallicity. The best evidence that TPs and the TDU occur in AGB stars is the formation 
of carbon st ars, which show C>0 at the ir surface as well as the presence of the radioactive 
element Tc (jWallerstein fc KnappI Il998l ) . This observed composition is explained by the 
combined effect of partial He burning during TPs (which produces more C than O) and 
the s process (which produces Tc), along with the TDU, which carries a fraction of the 
intershell to the stellar envelope. As mentioned above, HBB is present in AGB stars with 
masses higher than ~ 3 Mq at the metallicity considered here, converting C into N. 

The s process occurs in the intershell of AGB stars where He is abundant and (a,n) 
reactions can be efficier itly activated . For a detailed review on the s process in AGB stars 
we refer the reader to iBusso et al.l (119991 ). There are two possible neutron sources: the 
^^C{a,n)^^0 and the ^^Ne(a,n)^^Mg reactions. The ^^Ne neutron source is activated at 
temperatures higher than roughly 300xlO^K (hereafter MK). These high temperatures can 
only be reached inside the convective TPs, and mainly in AGB stars more massive than ~ 3 
Mq (jTruran &: Iben 1977 ). The ^^Ne is produced via double a captures on ^^N when the H- 
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burning ashes are engulfed in the TP. The resulting neutron flux occurs over a ti mescale of the 
order of a few years with relatively high neutro n densities, up to 10^^ n/cm^ (iGallino et al. 



19981 : iLugaro et al.l l2003l : Ivan Raai et al.l l201l[). The ^ ^C neutron source is instead acti- 



vated at temperatures as low as 90MK (jCameron 1955 ) but there is not enough ^^C in the 
H-burning ashes of AGB stars to make it an efficient neutron source. Without this low- 
temperature neutron source it is not possible to ex plain the s-pr o cess enhancements whic h 
have been observed in AGB stars since the 1950s (IMerrilll Il952l : ISmith fc LambertI Il986l ). 
This is because most of the observed s-process-enriched AGB stars are of relatively low mass 
(< 4 Mq) where the temperature is not high enough to activate the ^^Ne neutron source. 
Thus, it has become standard practice to artificially include the ^^C neutron source in AGB 
models of the s process. 

The most likely way to form the required ^^C is to allow some extra mixing of protons 
to occur following each TDU episode at the deepest extent of the convective envelope, where 
a sharp composition discontinuity forms at the interface between the H-rich envelope and 
the He-rich intershell. This extra m ixing may be physically ex plained by the effect of semi 



convection (jHoUowell fc Ibenlll988l). hydrodynamical overshoot (jHerwig]|2000l : ICristallo et al. 



2009bj), or gravity waves (jPenissenkov fc ToutI 120031 ) . However, there is no agreement yet 



on which of these proposed mechanisms is responsible for the extra mixing and thus we do 
not know what its detailed features are. The proton s are captured by the abundant ^^C, 



producing a layer rich in -^^C and ^^N (the ^^C pocket. iHoUowell &: Iben Il988t iGallino et al. 



19981 : iGorielv fc Mowlavill2000l : iLugaro et al.ll2003l : ICristallo et al.ll2009bh . The resulting 
usua lly burns under radiative conditions via "'^^C(a,n)^^0 before the onset of the following 



TP (iStraniero et al.lll995l : IGallino et al.lll998[ ). This burning releases a large number of free 
neutrons on a long timescale (around 10^ yr) and with low neutron densities < 10^ n/cm^. In 
the region of the pocket where ^^N is more abundant than ^^C, no s-process nucleosynthesis 
occurs owing to the dominance of the ^^N(n,p)^^C reaction over neutron captures by Fe seed 
nuclei and their progeny. In the intermediate-mass AGB stars that experience HBB, the 
formation of a ^^C pocket may be inhibited by proton captures occurring at the hot base of 
the c onvective envelope during the TDU, which produce ^^N and not ^^C (iGoriely &: Siess 
20041 ). Extremely deep TDU may also inhibit the activation of the ^^C(a,n) -'^^0 reaction b y 



penetrating into regions of the stellar core with a low abundance of "^He (iHerwigl l2004al ). 
While these processes do need further investigation, on first sight a dichotomy appears in 
models of the s process in AGB stars where the ^^C neutron source is activated in the low- 
mass models, while the ^^Ne neutron source is activated in the higher-mass models. At the 
metallicity we consi der here the mass at w hich the importance of the two neutron sources 
switches is ~ 3 Mq dGorielv fc Siessll2004f l. 



An interesting feature of the ^^C neutron source is that it is a primary neutron source 
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because it is produced from the H and He originally present in the stars. Since the total time- 
integrated neutron flux r from the ^'^C source is proportional to ^^C/Z, the 5 - process pre- 



dictio ns have a strong dependence on the metallicity of the star (iClaytonI 119881 : iBusso et al. 



pr 



200ll ). Specifically, with r increasing when decreasing Z , heavi er elements are produced at 



lower metallicities. This property allowed iGallino et al.l ( 1l998l ) to predict the existence of 
low- metallicity Pb-rich st ars, which was confirmed by observations of Pb-rich CEMP stars 



[e.g.. lVan Eck et al.ll2003l ). On the other hand, the ^^Ne neutron source is traditionally con- 
sidered to be a secondary neutron source, since it relies on the original CNO nuclei present 
in the star. However, at the low metallicities considered here ^^Ne also has an important 
primary component because the N abundance in the H-burning ashes increases significantly 
from its initial value. This is due to the effect of TDU carrying primary ^^C to the envelope, 
which is then converted into N by H burning. 

The ^^C and the ^^Ne neutron sources produce neutron fluxes and neutron-capture 
nucleosynthesis with different features. In summary, it has been shown that the ^^C source 
is responsible for the production of the bulk of the s-process elements in low-mas s AGB stars 
reaching, as mentioned above, up to Pb at low metallicities ( iGallino et al.lll998l ). Branching 
points on the s-process path are mostly closed during this neutron flux. The ^^Ne source 
instead produces smaller abundances of t he s-process elernents but it efficient ly activates 
branching points along the s-process path (lAbia et al.ll200ll : Ivan Raai et al.ll201ll ). Note that 
at low metallicity ^^Ne is both a neutron source and a neutr on poison via (n,7) reactions, 
leading to the produ ction of light elements such as Na and Mg (Herwig 2004bl : Karakas 2010 : 
Bisterzo et al.ll2010l ). Based on our models in § [3] we revise the different neutron-capture 



scenarios for AGB stars of low metallicity. 



Stellar Models 



The stellar evolutionary sequences were computed using the stars and t he Monash /Mou nt 



Stromlo (hereafter Stromlo) codes. The stars code was origina l ly written by 



Mixing is treated via a diffusion equation ( lEggletonl |1972| ) and convectiv e boundaries are 



EgdetonI fll97ll ) 



and has been updated by many authors fe.g. Pols et al. 1995 : Stancliffe fc Eldridgel 20091 ). 



determ ined via the Schwarzschild criterion. Th e opacities tables are from lEldridge fc Tout 



( I2OO4J ) and the nuclear reaction rates are from IStancliffe et al.l ( l2005l ). Models are evolved 
from the pre-main sequenc e to the AGB u sing 999 mesh points. Mass loss before the AGB 
phase is included using the iReimerj (I19751 ) presc ription, with rj = 0.4 . Mas s loss during the 
AGB phase is included using the prescription of IVassiliadis &: WoodI (119931 ). The computa- 



tion of the AGB phase with the stars code has been described in detail in IStancliffe et al. 
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( 120041 ) ■ We included an approximate contribution to the mole cular opacities (based on th e 
work of iMarigd |2002| ) that varies with envelope composition (jStancliffe fc Glebbeekll2008l ). 
Convective overshooting is not included at any stage in the evolution, and a mixing length 
of a = 2.0 is employed throughout. 

The Stromlo code uses the standard mixing-length theory for convective regions, with 
a mixing- length parameter a = 1.75, and determines the border of convective regions by 
applying the Schwarzschild criterion. Convective overshooting is not employed at any stage 
in the evolution, although the code sear c hes fo r a neutral border to the convective zone in 

This can have the effect of increasing the 



the same way as descr ibed in iLattanziol ( 



efficiency of the TDU (IFrost fc 



jattanzio 



1986 ) 



1996 



Karakas et al.ll2002l ). At high temperatures 



the OPAL opacities are used (llglesias fc Rogers! Il996l ). and at low temperatures (below 
T < lO^K) two different p rescriptions are employed. The models between 1 to 2.5Ma 



© 



(IKarakas fc Lattanzid 120071 ) include an approximate treatment fo r the molecul a r opa cities 
(in particul ar CN, CO, H2O, T iO) using the formulations from iBessell et al.l (119891 ) and 
corrected by lChiosi et al.l (119931 ). These fits do include some compositional dependence, but 
do not account for large variations in the envelope C/0 ratio or N abundance. In the model 
of 0.9 Mq and mo dels between 3 to 6 Mc :^ the code uses the (scaled-solar) lo w-temperature 
opacity tables from [Ferguson et al.l (120051 ) in place of t he iBessell et al.l ( 119891 ) fits. Mass loss 
before the AGE phase is included using the iReimersI (119751 ) prescription, with rj = 0.4. For 
the models bet ween 0.9 to 2.5 M^n mass loss during the AGB phase is included using the 



prescription of IVassiliadis &: Wood! ( 



9931) . For models of 3 to 6 Mq mass loss during the 



AGB is included using the iReimerj (119751 ) prescription with rj values that vary with mass 
(see discussion in lKarakasl[201o[ ). 



Two further Stromlo models of 2 Mq and 2.5 Mg^ of Z = 0.0001 were ca lculated using 



the low-temperature C and N-rich opacities from iLederer fc Aringei (2009). These new 
stellar evolutionary sequences were computed with IVassiliadis fc Wood! (119931 ) mass loss on 
the AGB and update d reaction rates for the ^^N(p, 7)^^0 reaction (the LUNA rate from 



Bemmerer et al.ll2006l ). and the NACRE rate for the triple-a process (lAngulo et al.lll999l ). 
These models were computed to quantitatively test the effect of the low-temperature opacity 
tables on the evolution of low-met allicity AGB stars, in particular to reconcile the number of 
TPs obtained between the Stromlo and the stars models. Many of the important features 
of the new models are the same as the old models, including the maximum depth of the 
TDU, however, the main change is a reduction in the number of TPs and hence the total 
amount of He-shell material dredged to the envelope over the AGB phase. One interesting 
feature of the new models is that they were both evolved to the white dwarf cooling track, 
where the final envelope masses are a few xlO~^ Mq, providing an excellent estimate of the 
total number of predicted TPs. 
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All our models were computed using a metallicity of [Fe/H]=— 2.3. The peak of the 
distribution of the metallicities of CEMP stars enhanced in neutron-capture ele ments is at 



[Fe/H ]=— 2.3 and the overall range cover —3.5 <[Fe/H]< —1. (see, e.g., Fig. 7 of lAoki et al. 



20071 ). Keeping this caveat in mind, we will compare our models to the composition of 



CEMP stars as a population, rather than to single objects. 

The main features of our stellar models are presented in Fig. [H and Tables [1] and |2l 
The tables include: the initial stellar mass, the total number of thermal pulses (TPs), which 
thermal pulses are followed by the TDU (TPs with TDU), the maximum temperature in the 
TPs followed by the TDU (Tgg^, also plotted in the bottom panel of Fig. [1]), the maximum 
TDU efficiency Amax, defined as the ratio of the mass scooped up to the envelope by the 
TDU over the mass growth of the H-exhausted core in the previous interpulse phase, the 
total mass dredged- up by TDU during the TP-AGB (total Mdred, also plotted in the upper 
panel of Fig. [1]), the mass of the envelope at the end of the computed evolution (final Menv), 
if HBB is at work or not, and the maximum temperature at the base of the convective 
envelope. The main features arising from the analysis of the two sets of models are: 

1. The number of TPs (Column 2) is smaller - roughly half at any given initial mass 
(Column 1) - in the STARS code than in the Stromlo code. This is because the stars 
code includes an approximation to the molecular opacities that varies with the com- 
position of the envelope. Th is is known to affect the mass-loss rates and shorten the 



stellar lifetime ( lMarigdl2002l ). The Stromlo 2 Mq and 2.5 Mq mod els computed using 



l ow-te mperature opacities with varied C and N abundances from iLederer fc Aringer 



( 120091 ) experienced ~ 40% fewer TPs thus confirming the effect of shortening the stel- 



lar lifetime when more realistic low-t emperature opacities a re used. The 2 Mq model 



of the same metallicity presented bv ICristallo et al.l (l2009bl ) uses the same molecular 



opacities from lLederer fc Aringerl (120091 ) . but a slightly different mass- loss law and has 
15 TPs, close to the number of TPs in our 2 Mq models. 

2. While the overall trend is for the number of TPs to increase with initial stellar mass 
(as a larger envelope means that it takes more time for the mass loss to erode it), 
models with masses lower than 1.25 Mq show instead an increase in the number of 
TPs. We obtain a higher total TDU mass and a higher intershell temperature in the 
0.9 Mq model than in the 1 Mq models. This result is most likely to be affected 
by uncertainties in the mass-loss rate, which are large, and needs to be investigated 
in more detail. In the lowest-mass models the initial pulses are too weak to trigger 
efficient TDU (Column 3), and thus C/0>1 is reached later during the evolution. 

3. The fact that the stars models experience fewer TPs with TDU (Column 3) than the 
Stromlo models leads to a smaller amount of dredged up material being mixed into the 
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envelope over the AGB phase (Columns 6). This is illustrated in Fig. [H The exception 
to this is the 1 Mq model. The total TDU mass does not change linearly with the 
number of TPs because it is also determined by the numerical treatment of the inner 
border of the convective envelope, which differs in the two codes. More efficient TDU 
is typically favored in the stars models (Column 5). 

4. The maximum temperature reached in TPs followed by the TDU (Column 4) is lower 
in the STARS models than in the Stromlo models (except for the STARS 1 Mq model, 
which has more TDU episodes). The ^^Ne(a,n)^^Mg reaction starts to operate at 
temperatures of ~300 MK, which is reached in the 1.9 Mq Stromlo model and in the 
2.5 Mq stars model. This quantity is also known to increase with the TP number. 

5. The final envelope mass is mostly determined by how long we managed to run the codes 
before encountering insurmountable convergence problems, usua lly during the final 



phase of very high mass loss that is typical of the prescription of IVassiliadis fc Wood 



( 119931 ) we used. It is usually easier to get closer to the end of the AGB phase with 
models of lower masses, although the intermediate-mass evolutionary sequences com- 
puted by the Stromlo code were al so evolve d to sm all (around 0.1 Mq) envelope masses 



probably due to the choice of t he iReimerd (119751 ) mass loss prescription on the AGB 



instead of IVassiliadis fc WoodI (119931 ). For these models we do not expect any further 
TPs and TDU episodes, except perhaps at the very tip of the AGB or during the 
post-AGB phase. 

6. The minimum mass at which we see the operation of HBB is 3 Mq, the same for both 
codes. HBB is very mild in this Stromlo model, with a maximum temperature at the 
base of the convective envelope ~40 MK (Column 8), while it is more efficient for the 
stars {T^^e —58 MK). HBB is already very efficient for the 3.5 Mq Stromlo model, 
with a maximum temperature ~76 MK. 

7. The Stromlo models with mass > 3 Mq reach temperatures in the He-burning shell 
up to 378 MK. These temperatures are high enough for an efficient activation of the 
^^Ne neutron source. The total TDU mass dredged up is similar in all intermediate- 
mass models, with a slight trend to increase with increasing stellar mass due to the 
increasing number of TPs. 



We note that most of the quantities presented in Tables H] and 12] are affected by uncer- 
tainties related to the treatment of the TDU and to the mass-loss rate. 
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2.1. Nucleosynthesis models 



To study the nucleosynthesis of the elements up to Bi we have used a post-processing 
code that takes as input the stellar evolutionary sequences described above. The stellar 
inputs are the temperature, density, and convective velocities (for convective regions) at 
each mass shell in the star as a function of time during the life of the star. Convective 
velocities are needed because the code includes both the changes due to nuclear reactions 
and those due to mixing in the equations to be solved for the abundances. An implicit 
method is used to solve N equations, N being the number of species included in the nuclear 



( 


2007 


), and 


Lugaro et al. 



Given that the number of species N determines how large the matrix to be solved is, 
in the interest of making the computation relatively fast (from a few days to a few weeks, 
depending on the initial stellar mass), we used for this work a basic s-process network, where 
we included a total of 320 species mostly along the valley of (3 stability, and 2,336 reactions. 
We based our nuclear reaction network on the JINA REACLIB database, as of May 2009. 
From this database, we updated the r ates for the neutron source r eactions: we took th e 
^^C{a,ny^O rate from lHeil et"aD fcoosh and the ^^lSle{a,nf^Mg from lKarakas et all J2006h . 



For the initial compos ition of most models we used the solar distribution of abundances 
from lAsplund et al.l (120091 ) scaled down to [Fe/H]=— 2.3. Solar ab undances of C, N, O, Ne, 
Mg, Si, S, Ar, and Fe are the pre-solar nebula values from Table 5 of lAsplund et al.l ( 120091 ): F 
is the meteoritic value of log e (F)© = 4.42 from Table 1 of the same paper (chosen because it 
has a lower uncertainty), and for many of the elements heavier than Fe we use the meteoritic 
values for the solar abundances (e.g., Sr, Eu, Pb). 

We also computed some models by varying the initial composition to: (1) the values 
predicted by galac tic chemical evolution m odels for solar neighborhood stars of metallicity 
[Fe/H]=— 2.3 from iKobayashi et al.l ( 120111 ). this case affects only the initial composition of 
the elements up to Zn; (2) the values obtained assuming that the star was born with an r- 
process enrichment and in some cases also an s-process enhancement. Enhancements in the 
initial r-process abundances may have arisen from the pollution from one or a few primordial 
SNell. Enhancements in the initial s-process abundances may have arisen from pollution 
from pri mordial massive stars producing s-process elements during core He and shell C 
burning (IPignatari et al.l 120081). or via the still unidentified "light element primary process" 
( iTravaglio et al.ll2004j : iMontes et al.ll2007l ). To compute these r and s-process enhancements 
we used the s and r contributions given in ISneden et al.l (120081 ) . 
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2.2. The inclusion of the ^^C pocket 



The methodology for the inclusion of the ^^C pocket deserves specific attention as the 
formation of this neutron source is still the main uncertainty in the s-process models. For 
reasons discussed above (§ ll.ip . we introduced a ^^C pocket for models of masses < 4 Mq, 
while we did not introduce it for the higher masses, except in one test case, where we 
introduced a ^'^C pocket in the Stromlo 5.5 Mq model. Also the 3 Mq model was computed 
with and without the inclusion of the ^^C neutron source to illustrate in detail the different 
effects of the two neutron sources. This mass lies where the two regimes switch: the s- 
process dominated by the ^^C neutron source and the s-process dominated by the ^^Ne 
neutron source. 

The inclusion of the ^^C neutron source was performed artificially during the post- 
processing by forcing the code to mix a small amount of protons from the envelope into the 
intershell at the deepest extent of each TDU episode. All the studies to date investigating the 
formation of the ^^C pocket by various mechanisms have found that the proton abundance 
in the intershell decreases monotonically. Hence, we apply the basic assumption that the 
proton abundance in the intershell decreases monotonically from the envelope value of ~ 0.7 
to a minimum value of 10"'' at a given point in mass located a t "Mmix" be low the base of 



the envelope. This method i s described in some more detail in iLugaro et al.l (120041 ) and is 
very similar to that used by iGoriely fc Mowlavil (120001 ) . We made the further choice that 
the proton abundance decreases exponentially, i.e., linearly in a logarithmic scale. Studies 
of the formation of the ^^C pocket have found profiles that can slightly differ from this 
basic assum ption, a. s well as from each other. For example, cor apare the profi l e show n 
in Fig. 4 of iHerwig (120001) to the profiles shown in Figs. 1-4 of ICristallo et al.l ( l2009bl ). 
Goriely Sz Mowlavi J2000l ) tested the effect of using a "slow" or a "fast" profile (see their 
Fig. 10) and found no major differences in the s-process distribution. 

We tested different values of our parameter Mmix, representing the e xtent in mass of the 
proton profile. The extent in mass of the result ing ^^C pocket is <Minix (IGoriely fc Mowlavi 
2000l : iLugaro et al.ll2003t ICristallo et al.ll2009al ). For some models of mass < 3 Mq we tested 



0, 0.0006 Mq, 0.001 Mq, 0.002 Mq, and 0.004 Mq. For the 3 Mq models we used 



'05 



'05 



0- 



'0 



0, 0.0005 Mq and 0.001 Mq. For the 5.5 Mq model, we tested M^i^=0 and 0.0001 



Mq. To keep the standard assumpt ion that the ^^C pocket should be a minor fraction 



(~1/10 - 1/20) of the intershell mass (IGallino et al.lll998l : iGorielv fc Mowlavil 120001 ) we had 

to decrease the value of M^ix with the stellar mass because the mass of the intershell is at 
least an order of magnitude smaller in the intermediate-mass AGB stars than in the lower- 
mass models: around 10~^ Mq in a 1.5 Mq model compared to around 10~^ Mq in a 6 Mp 
model. 



© 
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3. Results 

The results in terms of [X/Fe] for selected key elements (where X represent a given 
element) for all our models are presented in Tables [3] and [Hand Figs. [2] and [3d- The elemental 
abundances [X/Fe] in these tables are given at the stellar surface at the end of the computed 
evolution. The extended tables available on-line list, for each element X included in the 
network: log e(X), i.e., 12 + logio {N(X)/N(H)}, [X/H], [X/Fe], and [X/0], and X(i), which 
is the fraction of element i by mass at the surface. We also include the C/0 ratio and the 
isotopic ratios of C, N, O, and Mg. All elemental and isotopic ratios are number ratios. These 
quantities are tabulated at the stellar surface at the start of the evolution, after each TP, 
and at the end of the evolution, accounting for the radioactive decay of long-lived isotopes. 

We also present in the extended tables and in the figures of § H] the s-process in- 
dices light "Is" and heavy "hs". We choose the three main elements belonging to the 
first s-process peak Sr, Y, and Zr to define [ls/Fe] = ([Sr/Fe] + [Y/Fe] + [Zr/Fe])/3 and the 
three main elements belonging to the second s-process peak Ba, La, and Ce to define 
[hs/Fe] = ([Ba/Fe]-|-[La/Fe] + [Ce/Fe])/3. We choose these defin itions to compare our mod- 



els to the data from the compilation of iMasseron et al.l (120101 ). who report abundances for 
Ba, La, Ce, Eu, and PtB. We note that our definitions are different from those used by 



Bisterzo et al.l (120101 ) who took Y and Zr as the Is elements, and La, Nd, and Sm, as the 
hs elements. That choice was made in order to select the elements for which spectroscopic 
determinations are more reliabl^. In the exte nded on-line tables w e provide [Is/Fe] and 



[hs/Fe] as defined in this work and as defined by lBisterzo et al.l (l2010l ). for comparison. 



We have identified four different regimes of neutron captures in our models, which 
explain the different abundance distributions for the elements heavier than Fe in Fig. |51 
corresponding to the different [Ba/Sr] and [Pb/Ba] ratios listed in Table [3] and HI These 
regimes are described below, together with the typical stellar mass ranges at which they 
occur. The mass range in which each regime is at work is slightly different in the two sets 
of models, stars and Stromlo and in most models more than one regime is at work. In 
principle the final results depend on the interplay of the different regimes, though often a 
dominant situation can be identified. The regimes are: 



^We had numerical problems computing the nucleosynthesis for the STARS 1 AIq model, which we still 
have not manage to solve, so the results of this model are not included in the tables and figures. 

•^For data on the Is elements we searched the SAGA database ( Suda et al. 20081) . 

''Note that, according to Arlandini et al. ( 1999I) . about 80% of the solar abundances of Sr, Y, Zr, Ba, and 
Ce are produced by the s process, while La, Nd, and Sm have an important r-process contribution: 62%, 
56%, and 29% of their solar abundances are ascribed to the s process, respectively. 
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The ^'^Ne(a,np^Mg neutron source operates inside TPs when the temperature at the 
base of the TP reaches above 300 MK. From Fig. [T]it is clear that this neutron source 
operates when the initial stellar mass is higher than about 2 - 2.5 Mq. Two key 
features of this regime are that (i) the temperature is never high enough to completely 
burn the ^^Ne nuclei (at most the ^^Ne intershell abundance decreases by a factor of 
three) and (ii) the neutrons are distributed over the whole intershell by the convective 
TP. The resulting overall number of neutrons per Fe seed produces distributions where 
lighter s-process elements are favored with respect to heavier s-process elements. For 
this reason the [Ba/Sr] and [Pb/Ba] ratios are negative. The neutron density reaches 
up to ^ 10^^ n/cm'^, so branching points are affected but the overall distribution is 
a typical s-process distribution. In our models the activation of the ^^Ne neutron 
source in intermediate-m ass AGB stars does no t lead to an enhanced production of 



Eu, as was suggested by iMasseron et al.l (|2010[ ). The [Ba/Eu] ratio in these models 



remains close to unitjlf]. This is the value expected by the s process in equilibrium 
with a\NA = constant, where cr^ is the neutron-capture cross section of the isotope 
of mass A and A^a is its s-process abundance, considering that the neutron-capture 
cross section of, e.g., the neutron magic nucleus ^^^Ba is ~650 times lower than the 
neutron-capture cross section of, e.g., ^^^Eu, while its solar abundance is ~65 times 
highei0. 

2. The ^^C(a,n)^^0 neutron source operates in radiative conditions during the interpulse 
periods, with ^^C produced by the inclusion of a ^^C pocket (cases with Mmix 7^ 0). 
Typical models where this is the main regime of neutron captures are the models with 
initial masses between roughly 1.75 Mq and 3 Mq. Two key features of this regime are 
that (i) the ^^C nuclei burn completely in radiative conditions and (ii) the neutrons 
are captured locally in the thin ^'^C-pocket layer. The overall number of neutrons per 
Fe seed produces s-process distributions weighed towards Pb, with [Pb/Ba] ratios in 
some cases exceeding unity and 0< [Ba/Sr] <1. In the massive 5.5 Mq model the only 
effect of including the ^^C pocket is the increase in the production of Pb by a factor 
of ~50, so that this is the only massive AGB model where the final [Pb/Ba] ratio is 
positive. Efficient activation of the ^^C neutron source produces [Ba/Eu] ratios close 
to unity, which is typical of the s process and similar to what was found for Regime 



■^The [Ba/Eu] ratio decreases down to 0.11 only in the 3 Mq stars model computed without the inclusion 
of the ^^C pocket, but this is due to the fact that not much s process occurs here so that very little Ba and 
no Eu are produced in this model. 

^We can approximately apply the cta-^a = constant rule because the Eu abundance in s-process condi- 
tions follows that of Ba, the magic nuclei that precedes it in mass. 
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The ^^C(a,n)^^0 neutron source operates in convective conditions, with ^^C produced 
by the inclusion of a ^^C pocket (cases with M^aix 7^ 0). In this case the ^'^C in the 
^■^C pocket does not completely burn radiatively and is ingested in the following TP. 
This regime characterises the first two or three ^'^C pockets in the lowest mass models, 
e.g., in the 1.25 Mq, and 1.5 Mq models the first three ^^C pockets are engulfed in the 
following TPs while the abundance of ^^C is still higher than 10~^ (by number), and it 
is close to 10~^ for the first ^^C pocket. In the 1 Mq Stromlo model we find two TDU 
episodes and thus two ^^C pockets are included (cases with Mmix 7^ 0), both of which 
are ingested in the following TP while the abundance of ^^C is still higher than 10~^. 
The neutrons are released over material in the whole intershell region and the ^^C is 
mixed with the ^^N poison from the ^'^C pocket and from the H-burning ashes. Thus, 
the number of neutrons per Fe seed is lower with respect to the case when ^^C burns 
completely in radiative conditions. Consequently, the [Pb/Ba] ratio is also lower. 

The C(a,n)^^ O neutron source operates in convective conditions, with ^^C produced 
by the ingestion of a small number of protons in the TPs. This event occurs in the 
first few TPs in models with masses roughly less than 2.5 Mq. The convective TP 
reaches at most the point in mass where the H mass fraction is 0.01 and we never 
find a splitting of the convective TPs as is seen to occur at lower nietallici t ies and/or 



in other studies fe.g.. ICampbell fc Lattanzid l2008l : ICristallo et al.l l2009al : iLau et al. 



2OO9I : ISuda fc Fujimotd I2OIOI ). However, the modeling of the proton ingestion event 



strongly depends on the treatment of convection and the convective mixing scheme 
employed. In the Stromlo code we instantaneously mix convective regions at each 
iteration, which is not appropriate in cases such as proton ingestion episodes when 
mixing and bur ning timescales are comparable. In the stars code we use a diffusive 
mixing scheme ( lEggletonl Il972[ ). which would be able to resolve the splitting of the 
convective region if it occurred. On the other hand, the first 3D simulations of AGB 
proton ingestion episodes considered here suggest that an advective mixing scheme 
may be more appropriate than a diffusive mixing sch eme and that the splitting of the 
convective TP does not occur (jStancliffe et al.l 1201 ll ). With these caveats in mind, we 
find that proton ingestion episodes can be expected to lead to favorable conditions 
for s-process nucleosynthesis and to favour the production of the s-process elements 
lighter than Pb, for the same reasons related to Regime 3 discussed above. While for 
models with mass higher than ~ 1.9 Mq the effect of these events is dwarfed by the 
presence of the ^^C pocket, in the lower mass models the proton-ingestion events alone 
can result in significant production of s-process elements (cases with Minix=0). The 
results are relatively low [Pb/Ba] ratios and the usual s-process [Ba/Eu] ratio near 
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unity. The final [Ba/Sr] in these low-mass models depends on the detailed features of 
the proton-ingestion episodes, with ingestion of a smaller number of protons favouring 
the production of Sr, as well as on the amount of ^^C from the pocket that is left over 
to burn in the following TP (Regime 3). In some of the stars models the final result 
is a negative [Ba/Sr] ratio. Production of the light element Na due to the ^^Ne+p 
during the proton-ingestion events occurs in most of the low-mass models (cases with 
M„,i,=0). 

As outlined above, the Stromlo and stars models differ with regards to the masses at 
which the neutron-capture regimes are predominant. Thus, the relative distribution of the 
s-process elements (the [Ba/Sr] and [Pb/Ba] ratios) may be different for models of the same 
stellar mass (compare, e.g., the distributions plotted in Fig. [2] for the Stromlo and stars 
models of same mass 1.25 Mq.) The other main difference is that the absolute abundances 
with respect to Fe are typically smaller in the stars than in the Stromlo models. This is 
because, as discussed in § [21 the stars models have a smaller number of TDU episodes due 
to the inclusion of realistic molecular opacities which vary with the envelope composition 
and which result in a shorter AGB lifetime. The same effect, though to a smaller extent, 
appears when considering the Stromlo 2 Mq and 2.5 Mq models computed with the inclusion 
of molecular opacities. 

We changed the initial abundances in most of our nucleosynthesis post-processing models 
according to the different choices indicated in § 12.11 and did not find any remarkable difference 
in the any of the final abundances reported in Tables [3] and H] (as well as in Tabl^ discussed 
below). The only significant change occurred to [Eu/Fe] when we applied a significant r- 
process enhancement. This is obvious as Eu is a predominantly r-process element. For 
example, in the Stromlo 2 Mq model computed with Mmix=0.002 Mq, using an r-process 
contribution to the initial composition of [r/Fe]=0.0, 0.4, 1.0, and 2.0 dex, we obtained final 
[Eu/Fe]=1.51, 1.53, 1.60, and 2.08, respectively. In all cases, the [Ba/Fe] was unaffected, 
leading to final [Ba/Eu]=0.91, 0.89, 0.82, and 0.39, for the four cases respectively. 



3.1. The light elements: C, N, O, F, Ne, Na, and Mg 



The C, N, O, and F abundances in low-metallicity AGB stars h ave been alre a dy ex - 
tensively p resen ted by Herwiej (2004b), Karakas fc Lattanzio ( 2007 ). Lugaro et al. (2008), 



Stancliffd (|2009[ ) and iKarakad (120101 ) so here we only discuss them briefly. In Fig. E] we 



present the final [X/Fe] for the elements lighter than Fe for the same selected models of 
Fig. |2l In Table |5] we present the final [X/Fe] for C, N, O, F, and Ne for all our models and 
the indicated choices of Mmix- These light elements mostly do not change when varying this 
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free parameter: the only exceptions are for F and Ne. For example, F decreases by 0.41 dex 
and Ne by 0.45 dex in the Stromlo 1 Mq model when Minix=0 All our models are very rich 
in carbon and relatively poor in nitrogen, with typical [C/N] varying from ~ 1 dex, for the 
0.9 and 1 Mq models, to ^ 2 dex for models up to 3 Mq (Fig. [3]). From this stellar mass 
and above, HBB is activated and the nitrogen abundance increases and quickly overcomes 
the abundance of carbon. Oxygen is produced in our models, ranging from 0.42 to 1.34 dex. 
A typical product of AGB nucleosynthesis is fluorine, which we find to be strongly produced 
in our models: [F/Fe] is comparable to [C/Fe] up to the stellar mass where HBB is efficient. 
For higher masses fluorine is destroyed by proton captures. Neon is also a typical product of 
AGB star s, in the form of ^^Ne, and i ts production factor varies in our models from ~+l to 
~+3 dex (IKarakas fc Lattanzidl2003l ). This element is observable in planetary nebulae and 
can be used as a further probe of AGB nucleosynthesis. The Ne abundance was derived for 
the halo planetary nebula BoBn 1, which is also rich in C, N, an d F, and can be explained 
by low-metallicity AGB models of ~ 1.5 Mq (jOtsuka et al.ll2010l ). The production of F and 
Ne is primary in low-metallicity AGB stars because it i s driven by the TDU o f ^^C. Our 2 
Mq models produce very similar results to the model of I Crist alio et al.l (j2009bl ). 

The light elements Na and Mg require specific discussion because they depend on the 
abundance of ^^Ne, which in turn depends on the amount of TDU and therefore on the stellar 
mass, on the formation of the ^^C pocket, and on the overall neutron fiux. Our predictions 
for these elements are shown in Fig. [3] for selected models and listed in Tables [3] and |4] for all 
the models. Na is strongly produced in our models due to proton and neutron captures on 
^^Ne. For models without HBB, the increasing efficiency of the TDU with stellar mass means 
that the amount of primary ^^Ne and Na in the envelope also increase with mass up to ^ 2 
Mq\ this is illustrated in Fig. [31 In these low-mass models the ^^Na is sy nthesized by proton 
captures on ^^Ne during bot h the interpulse and convective TP (see also iGoriely fc Mowlavi 
2000l : ICristallo et al.ll2009br ). Neutron captures on 22Ne occur any time there are neutrons 



av ailable given the large abundance o f ^^Ne present in the intershell (see also the discussions 



m 



Karakasll2010l : iBisterzo et al.ll2010l ). In models of masses higher than around 2 M©, less 
Na is produced because the total amount of TDU as well as Mmix are smaller and because Na 
production by HBB is sensitive to mass. Na production requires higher temperatures than 
CNO cycling and takes place most efficiently at ^ 4 Mq; at higher masses the temperature 
at the base of the envelope is so high that the Na is effectively destroyed by proton captures. 
At the other end of our mass range, in the lowest-mass models 0.9 and 1 Mq, we also obtain a 
range of Na abundances, depending on the features of the proton ingestion and the inclusion 
of the ^^C pocket. The 0.9 Mq and 1.0 Mq Stromlo models are the only models to have 
[Na/Fe] significantly lower than [Ba/Fe]. In all the other models [Na/Fe] is comparable and 
even higher than [Ba/Fe]. 
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Magnesium is also significantly produced in our models when the mass is greater than ^ 
1.25 Mq. In the models with mass lower than ~ 2.5 Mq all the isotopes of Mg are produced 
via neutron captures in the ^^C pocket starting on ^^Na. For this reason the abundance of 
Mg follows the abundance of Na, keeping roughly to one order of magnitude lower than Na. 
In the models with mass higher than ~ 2.5 Mq, the heavier isotopes of Mg, ^^Mg and ^^Mg, 
are also produced by a-capture during the TP via the ^^Ne+a reactions. HBB in the most 
massive models leads to activation of the Mg-Al chain s, which produces ^^Mg and ^^Al and 



alters the composition of the He-shell prior to each TP (IKarakas fc Lattanzidl2003l ). In these 
higher mass models the abundance of Mg becomes comparable to, and even larger than, that 
of Na. 

In Fig. Hlwe show the Mg isotopic ratios at the stellar surface at the end of the computed 
evolution. These ratios vary by up to three orders of magnitude with the stellar mass. The 
derivation of the the Mg isotopic ratios from the spectra of CEMP stars would represent very 
strong constraints on the mass of the donor AGB star. We note that the models shown in 
Fig. Hlwere computed using the solar ratios of ^'^Mg/^^Mg=7.9 and ^^Mg/^^Mg=7.2. In se- 
lected models with mass > 1.5 Mq we tested th e effect of starting with ratios derived from the 



detailed Galactic chemical evolution models of iKobayashi et al.l ( 1201 ll ). i.e., ^'^Mg/^^Mg=158 
and 24Mg/26Mg=179. While this did not have any impact on the final ratios for the masses 
considered, the final compositions of stars with masses lower than 1.5 Mq may be more 
strongly effected by the initial composition. This is because these models experience less 
dredge-up and for this reason, the final Mg isotopic ratios are likely to be higher than given 
here. 



3.2. Comparison with previous models 



In Table [6] we p resen t a comparison of our results with those of ICristallo et al.l (l2009bf ) 
and iBisterzo et al.l (|2010| ) for the same element ratios as shown in Tables [3] and |H The 
results from the STARS and Cristallo's 2 Mq models show very good agreement, both in the 
light (see Table [5]) and in the heavy elements. Both these models were computed with the 
inclusion of molecular opacities. The Stromlo 2 Mq models both with and without inclusion 
of realistic molecular opacities produces higher abundances of both light and heavy elements 
in AGB models, with respect to Fe. Varying the mass of the ^^C pocket can also lead to 
significant changes in the predicted abundan ces of both light and heavy elements in AGB 
models. The model by lCristallo et al.l ( l2009bl ) includes the ^'^C pocket using a self-consistent 
approach based on convective overshooting motivated by hydrodynamical simulation of con- 
vective/radiative boundaries. These authors found that the extent in mass of the pocket 
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decreases with the TP number. In our models instead we kept Mmj ^ constant with the TP 



number. This explains why the abundances with respect to Fe by I Crist alio et al.l (l2009bf ) 
are closer to those computed with the smaller Mmix value in our models. In any case, the 
ratios [Ba/Sr], [Ba/Eu], and [Pb/Ba] are very similar in all three sets of models. This com- 
parison indicates that the s-process distributions obtained using our artificial treatment of 
the ^^C pocket are in agreement with those obtained using the self-consistent treatment of 
Cristallo et aP J2009bh . 



Bisterzo et al.l (120101 ) treated the neutron production in the ^^C pocket as a free param- 



eter and varied its efficiency with respect to their standard (ST) case, which reproduces the 
solar system s-process abundances at Z = 0.01. This approach results in a larger allowed 
range of variations of the resulting abundances. The results from th e Stromlo, STARS, an d 
Cristallo's models are typically included within the range predicted by lBisterzo et al.l ( 120101 ). 



Overall, Table [6] presents a relatively consistent set of results for neutron-capture nu- 
cleosynthesis in AGB stars of low metallicity. The main features of these models are the 
production of s-process elements up to +3 dex and a variety of s-process distributions, typi- 
cally with 0<[Ba/Sr]<1.0, [Ba/Eu]~0.9, and 0<[Pb/Ba]<2.0. [Ba/Eu] reaches lower values 
only if there is no t a large s-process enh ancement, i.e., [Ba/Fe]<1.0. The "radiative + con- 
vective" model of Goriely fc Siesj ( 2005 ) produces a heavy-element distribution in the range 
displayed in Table El with [Ba/Sr]~0.9, [Ba/Eu]~0.8, and [Pb/Ba]~1.0, see their Fig. 6. 
The p roduction of N a and Mg always increases with the stellar mass, something also seen 
in the iBisterzo et al.l ( 120101 ) models. Their 1.3 Mq models produce lower Na than the 1.25 
Mq models presented here because these models do not include production of Na via proton 
captures. 

Finally we find that, in spite of the fact that the mixing scheme strongly affects the 
proton-ingestion epi sodes, our results for the s process are qualitatively the same as the 
results presented by Cristallo et al. (j2009a ) for a 1.5 Mq star of metallicity half of that used 
in our models (Z = 5 x 1 0~^) which also experience a proton-ingestion episode. From Fig. 6 



of ICristallo et al.l (l2009al ) we see that after the first TDU episode, which records the effect 
of the proton-ingestion event, the s process favours the production of the lighter elements, 
just like our models. The final [Ba/Sr]~0.40, [Ba/Eu]^0.95, and [Pb/Ba]~0.85 are close to 
those obtained in our 1.25 Mq and 1.5 Mq models. 
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4. Comparison with the abundances of CEMP-s and CEMP-s/r stars 

4.1. [Ba/Fe] versus [Eu/Fe] 

Figure presents an overview of the observational data for CEMP stars and carbon- 
no rmal metal-poor stars of the rl and rll subclasses following the data and classification 



of iMasseron et al.l ( 120101 ) . In this and all the following figures we used different symbols 
for stars in different metallicity ranges, as indicated in the caption. There are no noticeable 
differences in the distribution of any of the plotted abundances when changing the metallicity, 
the only noticeable properties being that there are more CEMP-s than CEMP-s/r stars 
with [Fe/H]> -2, and more CEMP-no stars with [Fe/H]< -3 than in the other CEMP 
star classes. The different classes of CEMP stars appear to be fairly well distinguished in 
this plot. Most of the CEMP-no Eu data are upper limits, and the only CEMP-rll star, 
CS 22892-052, appears indistinguishable from the rest of the rll group in terms of its heavy 
element composition. 

The two groups of CEMP-s and CEMP-s/r appear to be distinct, with CEMP-s/r 
having higher [Eu/Fe] ratio in absolute terms, but also higher [Ba/Fe] on average. In Fig. [S] 
we plot the correlation lines of the two groups. The similar slope, around 1, indicates that 
for both classes the linear dependence would be preserved in a non-logarithmic scale. What 
counts here is the y-intercept, which corresponds to [Ba/Eu] and thus defines the slope of 
the linear correlation between the two variables in a non-logarithmic scale: in the case of 
CEMP-s stars the process responsible for their compositions produces on average of around 
400 times more Ba than Eu, with a spread of a factor of about 2; in the case of CEMP-s/r 
stars, instead, the process responsible for their composition produces on average of around 
200 times more Ba than Eu, with a spread of a factor of about 2. Not surprisingly, the r 
stars also sit on a correlation line with a slope of around 1 and y-intercept of around —0.6, 
corresponding to the solar-system r-process Ba/Eu ratio of about 10. In the case of CEMP- 
s the [Ba/Eu] value of +0.89 agrees very well with the typical [Ba/Eu] ratios reported in 
Tables [3l HJ and [6l while in the case of CEMP-s/r, the value of +0.6 does not agree with 
any AGB models producing [Ba/Fe] >1. 

The three CEMP-s/r with the lowest Ba abundance in their group may belong in- 
stead to the CEMP-s group. This seems likely for BS 16080-175, with [Ba/Fe] =1.55 and 
[Eu/Fe]=1.05, which can be reclassified as CEMP-s within the 0.07 dex uncertainties. The 
other two outliers BS 17436-058, with [Ba/Fe]=1.60 and [Eu/Fe]=1.17, and HD 187861, 
with [Ba/Fe]=1.39 and [Eu/Fe]=1.34, may be reclassified as CEMP-s within 2a of their 0.11 
and 0.20 dex [Eu/Fe] uncertainties. However, we note that [La/Fe]=1.73±0.2 in HD 187861 
indicates that it should be re-classified as a CEMP-s/r star. More observations and analysis 
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of these stars will help to clarify their status. 



Another interesting case is that of CS 30322-0 23 ([Ba/Fe]=0.52 and [Eu/Fe]=-0.63), 
which has been previously classified as a CEMP -s ( lAoki et al .1120071 ). and was included in 
the CEMP-low-s class bv iMasseron et all feoioh . This star has the highest [N/C] ~ +2 
and the lowest, and only negative, [Eu/Fe] ratio in the sample, as well as [Fe/H]^ —3.3, 
at the lower end of the CEMP-s range. Another three stars have [N/C] > 1: CS 22949- 
037, CS 22960-053, and CS 2 9528-041, which all ha ve [Fe/H] similar to CS 30322-023 and 
are classified as CEMP-no bv lMasseron et al.l fcoioh . However, CS 22960-053 and CS 29528- 
041 have [Ba/Fe]:^0 .9 but no Eu observa tions, so they may be belong to the C EMP-low- s 
class (as defined by IMasseron et al.ll2010l ). together with CS 30322-023. In fact. lAoki et e 
(2007) classified both CS 22960-053 and CS 30322-023 as CEMP-s stars following [pJan et al. 
fl2005h . 



Figure [6] includes a selection of our model predictions together with the data of Fig. |5l 
In this as in all the following figures the solid lines are a visual aid connecting the predicted 
envelope composition after each TDU episode. These prediction lines are not straight because 
the composition of the s-process material in the He intershell evolves with time and because 
the variables are shown on a logarithmic scale. The process of transferring and mixing 
AGB material onto the companion star would result in mixing lines connecting the AGB 
composition with the initial composition. These mixing lines closely follow the plotted AGB 
evolution lines and are not shown in the figure. We include in the figure the same stellar 
models plotted in Figs. |5] and [31 representing cases where the neutron-capture nucleosynthesis 
occurs in the different regimes described in § [31 All the plotted models start from a solar 
composition, [Ba/Fe] = [Eu/Fe] =0, except for the three 2 M© models from the Stromlo set 
with different initial composition: [r/Fe]=-|-0.4, [r/Fe]=-|-1.0, and [r/Fe]=-|-2.0 included in 
the lower panel of Fig. [SI 

All the models produce both Ba and Eu. The predictions lines follow the trend of the 
CEMP-s group. A small spread in the initial r-process composition up to [r/Fe]=±0.4 dex, 
which is withir i the composition obs erved in field stars of the corresponding metallicity (see, 
e.g.. Fig. 2 of iBisterzo et al.ll2011bl ). is enough to cover the data. CS 30322-023 indicates 
an initial negative [r/Fe]~ —0.6, which would also be within the observed range at the 
corresponding metallicity [Fe/H] =—3.3. On the other hand, AGB models do not produce 
the high [Eu/Fe] observed in CEMP-s/r stars. 

One solution to this problem may be to increase the initial [r/Fe] abundance by up 
to -|-2 dex under the assumption that the r-process enrichment is primordial in these sys- 
tems (bottom panel of Fig. This primordial enrichment can arise by assuming that the 
molecular clouds where CEMP-s/r stars formed were polluted by r-process elements from a 
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nearby SNII. This scenario provides a soluti on for the highest meas ured [Eu/Fe] abundances 
in CEMP-s/r stars as discussed in detail by lBisterzo et al.l (l2011b| ). Ho wever, it comes with 



three problems, som e of which have been discussed previously (see e.g. IJonsell et al.l 12006 
Lugaro et al.ll2009l ). (1) The initial [r/Fe] value does not affect the final [Ba/Fe] value 



which remains constant, hence the linear correlation observed between the Eu and Ba en- 
hancements in the CEMP-s/r sample is not reproduced. (2) The number of rll stars (around 
10) is smaller than the number of CEMP-s/r stars (about 30), so, it seems unlikely that 
CEMP-s/r stars evolved from rll stars in binary systems, because these should make up a 
fraction of all rll stars. This does not appear to be an observational bias since stellar surveys 
such as the Hamburg/ESO R-pro cess- Enhanced Stars (HERES) have specifically targeted rl 
and rll stars. It is also difficult to invoke supernova triggerin g as the cause for the forma tion 
of binaries: though early work indicated this link possible ( iVanhala &: Cameron! Il998l ). it 
is now understood that formation of binari es occurs v ery naturally whenever a turbulent 
velocity field is present in a molecular cloud ( lBatell2009l ). Furthermore the whole concept of 
super nova-triggered star formation is st ill controversial, both observationally and theoreti- 
cally ( iLeao et al.ll2009t lElmegreenI 120111 ). (3) The metallicity distribution is different for the 
two samples: [Fe/H] is centered at ^ —2.5 for CEMP-s/r stars and at ~ —2.8 for rll stars. 



4.2. [Is/hs] versus [Mg/hs] 

In this and the following subsections we compare our model predictions to the observa- 
tions using "intrinsic" indicators, i.e., elemental ratios that include only elements that are 
produced in AGB stars. Different from the [Ba/Fe] and [Eu/Fe] ratios considered above, 
where Fe is not produced in AGB models, the "intrinsic" ratios move away from their initial 
solar values towards their s-process values already after a small number of TPs (though this 
is less true for the 5 Mq model where the dilution in the envelope is much larger than in the 
other models). For the CEMP stars, this means that the intrinsic ratios are close to the AGB 
values even for relatively large dilutions of the AGB material onto the CEMP star. Hence, 
intrinsic ratios are, in a first approximation, independent of the TDU, the mass loss, the 
stellar lifetime, and the accretion and mixing processes on the binary companion. Instead, 
they mostly constrain the nucleosynthesis occurring in the deep layers of the star and in our 
context, the neutron source and the neutron flux for the s process. 

We also note that we consider here ratios of elements that are all significantly produced 
in AGB stars, so that variations in the initial composition of these elements due to the 
chemical evolution of the Galaxy do not have a significa nt effect on the AGB predictions. 



Using the initial composition from lKobayashi et al.l (1201 if ) (i.e., [F/Fe], [Na/Fe], and [Mg/Fe] 
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equal to —1.5, —0.2, and +0.3 dex, respectively) did not lead to modifications of the results 
plotted and discussed here in selected models with mass > 1.5 Mq. Only the [Mg/Fe] ratio of 
stars with masses lower than 1.5 Mq may be more strongly effected by the initial composition 
because these models experience less dredge-up. For these models, this effect would at most 
shift the [Mg/hs] ratios by +0.3 dex. 

Figure [7] shows the [Is/hs] and [Mg/hs] observed in CEMP-s and CEMP-s/r stars to- 
gether with our model predictions. The dispersion of the data is quite large, however, it 
clearly shows the presence of a linear correlation with slope ~ 1. The y- intercept ~ indi- 
cates that this correlation is most likely due to to variations in the amount of the hs elements 



Ba, La, and Ce. As n oted in previous studies (see, e.g.. Fig. 9 and Fig. 6 of lJonsell et al.ll2006 



Bisterzo et al.ll2011bl . respectively), the CEMP-s/r stars are separated from the CEMP-s in 
that they show [Is/hs] lower by at least 0.5 dex. Here we show that CEMP-s/r stars also 
show lower [Mg/hs] than CEMP-s. It seems obvious to ascribe these differences to the fact 
that CEMP-s/r stars have, on average, higher hs abundances than CEMP-s (Fig. |6]). 

All our models produce [Is/hs] > —1. Th is derives from the basic way in which the 



s process operates (see also IBusso et al.l 120011 ). As the neutron flux increases, first the Is 



elements, then the hs elements, and finally Pb are produced. As the flux goes through 
these peaks the abundances reach their equilibrium values, above which it is not possible to 
increase them. This is similar to the situation of buckets connected to each other. When 
water is poured into the flrst bucket, the bucket fllls up. Once it is full, the water starts 
moving into the second bucket and so on. The relative amount of water in the two buckets is 
only determined by their size (the neutron-capture cross sections in the case of the Is and hs 
elements) and not by how much more water is poured into the syste Hence, the minimum 
[Is/hs] ratio is the value expected when the s process is in equilibrium through the flrst and 
the second s-process peaks (with ctaNa = constant). This [Is/hs] ratio is roughly equal to —1 
because the neutron-capture cross section and the solar abundance of, e.g., ^^^Ba are about 
1.5 and 6 times lower (respectively) than those of, e.g., ^^Sr. On the other hand, the [Pb/hs] 
ratio cannot reach an equilibrium value because Pb is at the end of the s-process chain. The 
Pb abundance can continue to grow indeflnitely as more neutrons are made available. 

Unless we ascribe the overall mismatch to the large error bars, which still would not 
explain why CEMP-s/r stars are separated from CEMP-s stars in terms of both [Is/hs] and 
[Mg/hs], it appears that a fraction of CEMP-s/r stars do not carry a typical pure s-process 
signature in the relative abundances of their s-process elements. Together with Fig. |6l Fig. [7] 



^To finish with this simile, wc note that the third bucket, corresponding to the Pb peak in the s process, 
would have an infinite size. 
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demonstrates that (i) CEMP-s/r stars have the highest abundances of hs elements, and thus 
the lowest [Is/hs] and [Mg/hs] down to — 2|§, and (ii) our s-process models can match the 
observations well for CEMP-s stars, but cannot match those of CEMP-s/r stars, not only 
in terms of their Eu abundance, but also in terms of their Is and hs abundances. 

Focusing on the CEMP-s stars only, they appear to be overall better matched by the 
models with the ^^C neutron source burning radiatively. The 5 Mq model, which could not 
be ruled out as a match to the observations on the basis of the Ba and Eu enhancements 
(Fig. E]), produces too much Mg and Is with respect to hs to be able to match any of the 
CEMP-s and CEMP-s/r data. Inspection of the results for the 5.5 Mq model in Table H] 
shows that inclusion of the ^^C pocket does not solve this problem. On the other hand, 
the composition of CS 30322-023, the only CEMP-s with [Mg/hs]>0, is compatible with the 
models of 3 Mq and 5.5 Mq computed with the inclusion of the ^'^C pocket. 



4.3. [Pb/hs] versus [Mg/hs] 

Figure [8] presents the [Pb/hs] and [Mg/hs] observed in CEMP-s and CEMP-s/r stars 
together with our model predictions. The dispersion of the data is large and there is no ob- 
vious correlation between the two ratios. While the two classes are different in all the other 
abundance ratios considered so far, the distribution of [Pb/hs] in CEMP-s and CEMP- 
s/r stars covers the same range. The models with the ^^C neutron source burning radia- 
tively provide a good match to CEMP-s stars, where the spread in the [Pb/hs] ratio can 
be explained by varying the T DU number or, equivalently, the initial stellar mass (see also 



Bonacic Marinovic et al.l 120071 ). The models where the ^^C neutron source burns convec- 
tively can provide a match for the lower end of the [Pb/hs] ratios observed in CEMP-s/r 
stars, however, as shown in the previous section, they do not match the [Is/hs] values of 
CEMP-s/r. Even though the error bars are large, CEMP-s/r stars with high [Pb/hs] ratios 
are not matched by any models. 

We note that the 5 Mq model does not cover any of the data points. As in the case of the 
[Is/hs] versus [Mg/hs] the composition of CS 30322-023 can be matched by intermediate-mass 
AGB models (e.g., 3 Mq and 5.5 Mq) if a ^'^C pocket is included. 



8 Specifically, the CEMP-s/r star with [ls/hs]~ -2 is HE 0212-0557 observed bv lCohen et al.l (|2006l ). The 
hs value of this stars is based on values of [Ba/Fe], [Ce/Fe], and [La/Fe] all very similar to each other, while 
Is is based on [Y/Fe]=0.55, derived from three lines, and [Sr/Fe]=— 0.05, derived from one line. Removing 
Sr form the calculation of Is still results in a very low [ls/hs]=— 1.69. 
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4.4. [Is/hs] and [Pb/hs] versus [Na/hs] and [F/hs] 

Figures [9] and [TO] present the [Na/hs] and [F/hs] ratios plotted versus the intrinsic s- 
process indicators, [Is/hs] and [Pb/hs]. There are not enough data to be able to determine 
if CEMP-s and CEMP-s/r stars are distinguished in their [Na/hs] and [F/hs] abundance 
ratios. The models where the ^^C neutron source burns radiatively, which provide a good 
match to the overall composition of CEMP-s stars in terms of their [Mg/hs], [Is/hs], and 
[Pb/hs] compositions, typically produce too much Na and F with respect to the hs elements 
to match any of the observational data. One exception is the case of CS 30322-023, the only 
CEMP-s with [Na/hs] >0: this star is again distinguished from the rest of the CEMP-s stars 
and as mentioned above may carry the signature of a more massive AGB star. A better 
overall match to [Na/hs] and [F/hs] is provided by the low-mass Stromlo models where the 
^^C neutron source burns convectively, however, these models do not cover the corresponding 
[Pb/hs] (and [Is/hs]) ratios. 

As discussed in § 13.11 the predicted abundance of Na is made by both neutron and 
proton captures on ^^Ne. The fraction made by neutron captures is linked to the production 
of Mg, and the models where the ^^C neutron source burns radiatively provide a good 
match to the [Mg/hs] observed in CEMP-s. The fraction made via proton captures, on the 
other hand, is disconnected from the production of Mg, but depends on the details of the 
proton-ingestion episodes, the profile of protons introduced to make the ^^C pocket, and th e 



proton-capture reaction rates of ^^Ne and ^^Na (which are uncertain, see llliadis et al.ll2010l ). 
All these uncertainties need to be tested to see if it is possible to match the observed [Na/hs] 
and [Mg/hs] ratios, which indicate that the production of Na via proton captures should 
be minimised. During the proton diffusion leading to the formation of the ^^C pocket, a 
pr oton profile weighed tow ards a low proton abundance (e.g., the "fast" profile in Fig. 10 



of iGoriely fc Mowlavil 120001 ) may contribute to the solution of this problem. This possibility 
will be investigated in detail in a forthcoming study and has the potential to provide a strong 
constraint on the formation of the ^^C pocket. A lower value of the ^^Ne(p,7)^^Na reaction 
rate or a higher value of the ^^Na(p,a)^°Ne reaction rate could also help lower the predicted 
Na abundances. 

The predicted [F/hs] ratios are also higher than the observed upper limits. The operation 
of some form of deep mixing of the envelope material of the AGB companion star down 
to regions where H burning occurs, which is also invoked to explain the observed high N 



enhancements and low ^^C/^^C ratios (see, e.g.. IStancliffell2010l ) may contribute to decreasing 
the F abundance in these stars. Further investigations are required to test this idea. Note 
that in the 5 model with HBB, F is strongly depleted by proton captures. 
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5. Summary and conclusions 

We have presented a large set of stellar models at metallicity [Fe/H] =— 2.3 from two 
different evolutionary codes for which we have computed detailed s-process nucleosynthesis. 
For the first time we have presented abundance predictions for the elements heavier than iron 
from stars with masses lower than 1.3 Mq and above 2 Mq at this metallicity. The most 
up-to-date models including the effect of molecular opacities that vary with the envelope 
composition during the AGB evolution experience less TPs, less TDU episodes, and produce 
overall lower enhancements of the s-process elements with respect to Fe. 

We have found that there are four regimes of neutron-capture processes and the domi- 
nant regime depends on the initial stellar mass: for the highest mass models the ^^Ne neutron 
source is dominant and the inclusion of a ^^C pocket where ^^C burns radiatively has an effect 
only on the production of Pb; for models of mass < 3.5 M©, ^^C is the main neutron source. 
For masses ~ 2 Mq, ^^C burns radiatively in the ^^C pocket during the interpulse periods. 
For lower masses, ^^C also burns convectively as the first few ^^C pocket are engulfed in the 
following TPs before ^^C is completely burnt in the interpulse period. For the lowest masses, 
proton-ingestion episodes associated with the first few TPs also produce ^'^C, which burns 
convectively inside the TPs. Neutrons released during convective TPs produce lower [Pb/hs] 
ratios than neutrons released in the ^'^C pocket under radiative conditions, and higher [Is/hs] 
ratios if a smaller number of protons are ingested. It should be kept in mind that proton 
ingestion episodes are extremely sensitive to the convective mixing scheme adopted. Our 
models provide some insight into the effect of thes e episodes on the 5-pro cess distribution. 



and are in agreement with the results prese nted bv ICristallo et a. 



detailed 3D models of the kind presented by IStancliffe et al.l (1201 if ) are needed to accurately 



(l2009al ). However, more 



ascertain the effect of proton-ingestion episodes on the s-process in AGB stars. The pro- 
duction of Mg and Na increases with the stellar mass as there is more primary ^^Ne seed 
available. 



From comparison with the observational data-set of CEMP stars compiled by lMasseron et al 



(120101 ) we have derived that our AGB models, especially those with radiative ^'^C burning, 
present a very good match to the composition of CEMP-s stars. On the other hand, they 
cannot match the compositions of CEMP-s/r stars, since they cannot produce the observed 
Eu enhancements up to one order of magnitude higher than in CEMP-s stars. If the s 
process and r process components, making Ba and Eu, respectively, in CEMP-s/r stars are 
completely independent, as in the scenario where the r enrichment is primordial, it is diffi- 
cult to explain why the observed Ba and Eu are correlated. This is because for any given 
initial [r/Fe] abundance pattern imposed, the AGB models produce the same final [Ba/Fe] 
abundance. Since this conclusion relies on the observed correlation between Ba and Eu, 
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it is urgent to confirm and establish this correlation using a larger number of CEMP-s/r 
star abundances derived from high-resolution spectroscopy and for a larger number of hs 
elements. 

Not only do CEMP-s/r and CEMP-s stars differ in that CEMP-s/r stars have more Eu 
than CEMP-s stars in absolute terms, but they also differ in that CEMP-s/r stars have on 
average more Ba (and other hs elements) than CEMP-s stars. This behaviour is not readily 
explained if the s-process component in CEMP-s and CEMP-s/r stars come from the same 
AGB star sources. Furthermore, the higher Ba (hs) abundances are not accompanied by 
higher abundances in Mg and in the Is elements, as it is expected in our models. In fact, our 
AGB models always produce [ls/hs]> —1, which is too high to match the distribution of the 
s-process elements in CEMP-s/r with [Is/hs] down to ~ —2. On the other hand the higher 
Ba abundances are accompanied by higher Pb abundances, as expected in the s process. In 
summary, the composition of CEMP-s/r stars is characterised by large enhancements of the 
hs elements and of Eu (correlated with Ba: ~ 1/200 Ba), as well as of Pb, but not of the Is 
elements. It remains to be seen if these compositions could be achieved by a "s/r" neutron- 
capture process, which may be a single neutron-capture process with features in-between or 
a superposition of the s and r processes. Parametric models of neutron-capture processes 
will be a first step in establishing if this is a possibility. 

If the Ba versus Eu correlation in CEMP-s/r stars is confirmed and if it is not due to 
a new "s/r" neutron-capture process, but to the s process and the r processes occurring 
in AGB stars and SNIIe, respectively, then these sites must be somehow correlated to each 
other. A probably unlikely scenario, which we provide here as an illustrative example, is that 
of a stable triple stellar system comprising the CEMP-s/r star observed today, a former AGB 
star, and a former SNII. In this case the orbital parameters such as the distance between 
the stars, which affect the mass transfer, are correlated, probably resulting in correlated Ba 
and Eu abundances. 

While the AGB models produce a good match to the observed [Mg/hs] ratios in CEMP- 
s stars, they typically produce too high [Na/hs] and [F/hs] ratios. This may depend on the 
choice of the proton profile leading to the formation of the ^^C pocket, the proton-capture 
reaction rates involved, and the occurrence of extra-mixing processes. Investigation of these 
issues may result in important constraints on these uncertain stellar and nuclear inputs. 

Finally, CS 30322-023 may carry the signature of nucleosynthesi s in intermediat e -mass 



AGB stars. This would be difficult to reconcile with the suggestion of iMasseron et al.l (120061 ) 
that this star is an intrinsic AGB object based on the low \ogg = —0.3. This would make 
this star an extremely peculiar object born less than a half a Gyr ago, but with an extremely 
low metallicity typical of the early Universe. Our conclusions are more in line with the 
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higher value of log (7 = 1 reported by lAoki et al.l (120071 ). not indicative of an intrinsic AGB 
star. From the preliminary comparison made here the composition of CS 30322-023 indicates 
that both HBB and the ^^C pocket were operating in low-metallicity massive AGB stars. 
However, before drawing conclusions we need to compare the composition of this star with 
predictions from a model of the appropriate metallicity of [Fe/H] =— 3.3. Observations of 
Eu in CS 22960-053 and CS 29528-041 may reveal that these two stars are also CEMP-low-s 
with [N/C]>1 creating a new class together with CS 30322-023. Future studies will focus on 
comparison of single stellar objects with our AGB model predictions. 
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Mass (Mq) 



Fig. 1. — The total mass dredged- up by the TDU (top panel) and the maximum temperature 
at the base of TPs (bottom panel) as a function of the initial stellar mass for our AGB models. 
The change in the trend for the Stromlo models between 2.5 Mq and 3 Mq i s caused by the 
change in the adopted mass- loss law (see the text for details). The use of the iReimerd (119751 ) 
mass-loss prescription for stars with masses 3 Mq - 6 Mq has the effect of shortening the 
AGB lifetime and reducing the number of TDU episodes. 
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Fig. 2. — The [X/Fe] ratio as function of the atomic number Z for elements heavier than Fe 
for selected Stromlo (upper panel) and STARS (lower panel). Data are surface abundances 
after the last computed TP. The selected models represent cases where the neutron-capture 
nucleosynthesis occurs according to the four regimes: The 5 Mq model represents the case 
where all neutrons are released by the ^^Ne neutron source in the TPs (Regime 1). The 2 
Mq models represent the regime where most neutrons are released under radiative conditions 
during the interpulse periods by the -"^'^C neutron source (Regime 2). The STARS 1.25 AIq 
model represent a case where neutrons released in radiative conditions by the ^'^C neutron 
source primarily determine [Pb/Ba], while [Ba/Sr] is also affected by ^^C burning convectively 
(Regimes 2, 3, and 4). The Stromlo 1.0 Mq and 1.25 Mq models represent the regime where 
most of the neutrons are released in the convective TPs by the ^^C neutron source (Regimes 
3 and 4). All the other models show behaviour in between the selected models and are not 
included for the sake of clarity (all the data can be found in the on-line tables). 
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Fig. 3. — Same as Fig. [2] but showing the [X/Fe] ratios as a function of the atomic number 
Z for elements hghter than Fe. Stromlo models are shown in the upper panel and stars 
models in the lower panel. 
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Fig. 4. — Mg isotopic ratios (by number) as function of the initial stellar mass: open symbols 
represent Stromlo models, black symbols represent STARS models. Squares are for models 
with Minix=0 and triangles are for models with Minix=0.002 Mq. 
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Fig. 5. — Overview of the observational data for CEMP stars (colour-coded as in the legend) 
an d rl (black crosses) ari d rll (open star symbols) stars following the data and classification 
of iMasseron et al.l (l2010l ) except that, for sake of simplicity, we have included the CEMP- 
low-s in the CEMP-s group. For CEMP stars, triangles represent stars with [Fe/H] > —2, 
squares represent stars with —3 <[Fe/H]< —2, and circles represent stars with [Fe/H]< —3. 
The correlation lines going through CEMP-s (short-dashed line corresponding to [Ba/Fe] = 
0.98 [Eu/Fe] -|- 0.89, with a correlation coefficient r=0.58) and CEMP-s/r (long-dashed line 
corresponding to [Ba/Fe] = 0.92 [Eu/Fe]+ 0.60, with correlation coefficient r=0.59) are also 
plotted. (A color version of this figure is available in the online journal.) 
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[Eu/Fe] 



Fig. 6. — A selection of our model predictions is plotted together with the CEMP data of 
Fig. |5l The predictions lines from our models represent the evolution of the stellar surface 
composition of the primary AGB stars and each symbol on the lines represents a TP. Top 
panel: open symbols refer to the Stromlo models and black symbols refer to the stars 
models. Stellar masses are indicated in the legend. Bottom panel: open symbols refer to 
Stromlo models of 2 Mq and different initial compositions as indicated in the legend. (A 
color version of this figure is available in the online journal.) 
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Fig. 7. — Same as Fig. [6] but for [Is /hs] versu s [Mg/hs]. For data on the Is elements we 
searched the SAGA database ( Suda et al. 20081 ). (A color version of this figure is available 
in the online journal.) 
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Fig. 8. — Same as Fig. |6]but for [Pb/hs] versus [Mg/hs]. (A color version of this figure is 
available in the online journal.) 
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Fig. 9. — Same as Fig.|6]but for [Is/hs] versus [Na/hs] ( t op pa nel) and [Pb/hs] versus [Na/hs] 
(bottom panel). Data for Na are from iLucatello et al.l ( 120111 ). (A color version of this figure 
is available in the online journal.) 



-43 - 



1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




CEMP-s 
CEMP-s/r 




I I I I I I I I I I I I I I I I I L 

-2-10 1 

[F/hs] 



Fig. 10. — Same as Fig. |6]but for [Is/ hs] versus [F/hs] (top p anel) and [Pb/hs] versus [F/hs] 
(bottom panel). Data for F are from iLucatello et al.l ( 201ll ). (A color version of this figure 
is available in the online journal.) 



Tabic 1: Details of stellar models computed with the STARS code. 



Mass (Mq) 


TPs 


TPs with TDU 


TJ^I^ (MK) 




total Mdred {Mq) 


final Menv (Mq) 


HBB 


T:^Z (MK) 


1.0 


13 


3-8 


273 


0.127 


0.006 


0.016 


No 


1.64 


1.25 


9 


2-7 


267 


0.431 


0.021 


0.015 


No 


2.75 


1.5 


9 


2-8 


268 


0.676 


0.040 


0.408 


No 


3.97 


1.75 


11 


2-10 


276 


0.846 


0.074 


0.08 


No 


6.05 


2.0 


13 


2-12 


283 


0.971 


0.098 


0.468 


No 


9.79 


2.25 


13 


2-13 


294 


1.060 


0.082 


0.584 


No 


22.04 


2.5 


16 


2-16 


298 


1.083 


0.107 


0.789 


No 


32.05 


3.0 


16 


2-16 


331 


1.116 


0.047 


0.979 


Yes 


56.76 



Table 2: Details of stellar models computed with the Stromlo 


code. 








Mass (Mq) 


TPs 


TPs with TDU 


Tg,^^ (MK) 


Amax total Mdred (Mq) 


final Menv (Mq) 


HBB 


^max (MK) 


0.9 


38 


4-20 


278'^ 


0.14 


0.012 


0.029 


No 


1.8 


1.0 


26 


5-6 


261^ 


0.18 


0.003 


0.0085 


No 


1.5 


1.25 


16 


3-13 


274^ 


0.25 


0.020 


0.018 


No 


2.3 


1.5 


18 


3-17 


282 


0.51 


0.059 


0.022 


No 


3.4 


1.75 


20 


3-20 


294 


0.71 


0.132 


0.029 


No 


5.4 


1.9 


24 


3-23 


300 


0.80 


0.169 


0.029 


No 


7.4 


2.0 


26 


2-26 


307 


0.85 


0.218 


0.038 


No 


9.0 




19'^ 


2-18 


305 


0.85 


0.133 


0.0^ 


No 


6.6 


2.25 


33 


3-33 


320 


0.94 


0.292 


0.162 


No 


15.7 


2.5 


30 


2-30 


351 


1.01 


0.283 


0.413 


No 


28.4 




21'' 


1-20 


341 


0.99 


0.165 


0.0^ 


No 


16.3 


3.0 


20 


2-20 


336 


1.02 


0.106 


0.119 


Yes 


39.6 


3.5 


27 


2-27 


348 


1.04 


0.103 


0.136 


Yes 


75.9 


4.0 


37 


2-37 


352 


0.97 


0.128 


0.133 


Yes 


82.9 


4.5 


41 


2-41 


352 


0.96 


0.110 


0.085 


Yes 


87.2 


5.0 


56 


2-56 


363 


0.96 


0.106 


0.124 


Yes 


92.5 


5.5 


77 


2-77 


369 


0.94 


0.124 


0.085 


Yes 


98.2 


6.0 


109 


2-109 


378 


0.93 


0.127 


0.183 


Yes 


105 



"This temperature refers to the TP followed by TDU. The absolute maximum temperature was 302 MK. ^This 
temperature refers to the TP followed by TDU. The absolute maximum temperature was 298 MK. ^This temperature 
refers to the TP followed by TDU. The ab solute maximum te r npera ture was 277 MK. '^This model was computed 
using the C and N-rich opacity tables from iLederer fc Aringed (120091 ). ^This model was evolved to the white dwarf 
cooling track. The final envelope mass was a few 10~^ Mq, which is the remaining H-shell. 



Table 3: Selected key elemental abundances [X/Fe] at the stellar surface at the end of the computed evolution for all 
our STARS models with a scaled solar initial composition from Asplund et al (2009). 



Mass (Mq) 


Mmix (Mo) 


[Na/Fe] 

L / J 


[Mg/Fe 

L O/ 


Sr/Fe 


[Ba/Fe 

L / 


Eu/Fe] 

/ J 


[Pb/Fe] 

L / J 


[Ba/Sr] 

L / J 


[Ba/Eu] 

L / J 


[Pb/Ba] 

L / J 


1.25 





0.65 


0.08 


0.36 


0.02 


-0.01 


0.02 


-0.34 


0.03 


0.00 




0.002 


1.62 


0.51 


1.65 


1.59 


0.69 


2.75 


-0.06 


0.90 


1.16 


1.5 





0.92 


0.21 


0.07 


0.02 


-0.01 


0.01 


-0.05 


0.03 


-0.01 




0.0006 


1.31 


0.42 


1.46 


1.38 


0.47 


2.49 


-0.08 


0.91 


1.11 




0.002 


1.86 


0.76 


1.64 


1.78 


0.85 


2.91 


0.14 


0.93 


1.13 




004 


2 28 


1 03 


1 34 


2 03 


1 1 1 


3 1 


69 


92 

yj .u £j 


1 07 

^..\J \ 


1.75 





1.58 


0.68 


0.10 


0.04 


-0.02 


0.02 


-0.06 


0.12 


-0.02 




0.002 


2.24 


1.15 


1.25 


1.94 


1.03 


3.05 


0.69 


0.91 


1.11 


2.0 





1.65 


0.92 


0.30 


0.64 


0.08 


0.44 


0.34 


0.56 


-0.20 




0.0006 


1.90 


1.09 


0.90 


1.52 


0.67 


2.69 


0.62 


0.85 


1.17 




0.002 


2.26 


1.32 


1.34 


2.04 


1.12 


3.03 


0.70 


0.92 


0.99 




0.004 


2.53 


1.52 


1.59 


2.31 


1.38 


3.12 


0.72 


0.93 


0.81 


2.25 


0.002 


1.86 


1.38 


1.18 


1.85 


0.92 


2.72 


0.67 


0.93 


0.87 


2.5 


0.001 


1.40 


1.24 


0.85 


1.47 


0.60 


2.48 


0.62 


0.87 


1.01 




0.002 


1.64 


1.42 


1.06 


1.80 


0.88 


2.59 


0.74 


0.92 


0.79 


3.0 





0.67 


0.69 


0.47 


0.11 


0. 


0.02 


-0.36 


0.11 


-0.09 




0.0005 


1.23 


1.19 


0.83 


1.25 


0.42 


2.24 


0.42 


0.83 


0.99 
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This preprint was prepared with the A AS macros v5.2. 



Table 4: Selected key elemental abundances [X/Fe] at the stellar surface at the end of the computed evolution for all 
our Stromlo models with a scaled solar initial composition from Asplund et al. (2009). 



Mass (Mq) 


M„,ix (Mq) 


[Na/Fe] 


[Mg/Fe] 


;Sr/Fe] 


[Ba/Fe] 


[Eu/Fe] 


[Pb/Fe] 


[Ba/Sr] 


[Ba/Eu] 


[Pb/Ba] 


0.9 





1.27 


0.13 


1.55 


2.24 


1.05 


1.85 


0.69 


1.19 


-0.39 




0.002 


2.07 


0.67 


1.56 


2.35 


1.31 


3.02 


0.79 


1.04 


0.67 


1.0 





0.30 


0.02 


1.14 


2.05 


1.08 


2.29 


0.91 


0.97 


0.24 




0.002 


0.64 


0.02 


1.10 


2.03 


1.11 


2.46 


0.93 


0.92 


0.43 


1.25 





0.60 


0.05 


1.37 


0.76 


0.03 


0.07 


-0.61 


0.73 


-0.69 




0.002 


1.68 


0.50 


1.51 


2.17 


1.08 


2.76 


0.66 


1.09 


0.59 


1.5 





1.51 


0.31 


1.57 


2.12 


0.97 


1.80 


0.55 


1.15 


-0.32 




0.0006 


1.78 


0.55 


1.59 


2.38 


1.37 


2.88 


0.79 


1.01 


0.50 




0.002 


2.30 


0.91 


1.64 


2.45 


1.47 


3.15 


0.81 


0.98 


0.70 




0.004 


2.68 


1.20 


1.78 


2.51 


1.53 


3.20 


0.73 


0.98 


0.69 


1.75 





2.30 


0.85 


1.83 


2.27 


0.96 


1.55 


0.44 


1.31 


-0.72 




0.002 


2.75 


1.28 


1.93 


2.54 


1.46 


3.22 


0.61 


1.08 


0.68 


1.9 


0.002 


2.90 


1.45 


1.90 


2.42 


1.45 


3.23 


0.52 


0.97 


0.81 


2.0 





2.73 


1.35 


0.62 


0.12 


-0.05 


0.04 


-0.50 


0.17 


-0.08 






2.09 


0.91 


0.57 


0.16 


-0.03 


0.03 


-0.41 


0.19 


-0.13 




0.0006 


2.81 


1.45 


1.51 


1.89 


0.98 


2.92 


0.38 


0.91 


1.03 




0.002 


3.02 


1.62 


1.78 


2.42 


1.51 


3.24 


0.64 


0.91 


0.82 




0.002^^ 


2.57 


1.30 


1.56 


2.21 


1.30 


3.17 


0.65 


0.91 


0.96 




0.004 


3.21 


1.80 


1.98 


2.65 


1.68 


3.28 


0.67 


0.97 


0.63 


2.25 


0.002 


3.07 


1.94 


1.81 


2.51 


1.55 


3.20 


0.70 


0.96 


0.69 


2.5 


0.001 


2.60 


2.05 


1.47 


2.11 


1.18 


3.04 


0.64 


0.93 


0.93 




0.002 


2.71 


2.13 


1.72 


2.40 


1.44 


3.14 


0.68 


0.96 


0.74 




0.002^^ 


2.28 


1.67 


1.50 


2.21 


1.26 


3.04 


0.71 


0.95 


0.83 



"Model computed with low-temperature opacities from lLederer fc Aringeri ( l2009l ). 



Table 4: Continues. 



Mass (Mq) 


M„,ix (Mq) 


[Na/Fe] 


[Mg/Fe] 


[Sr/Fe] 


[Ba/Fe] 


[Eu/Fe] 


[Pb/Fe] 


[Ba/Sr] 


[Ba/Eu] 


[Pb/Ba] 


3.0 





1.91 


1.72 


1.40 


0.79 


0.04 


0.19 


-0.61 


0.75 


-0.60 




0.0005 


2.08 


1.87 


1.43 


1.76 


0.78 


2.77 


0.33 


0.98 


1.01 




0.001 


2.17 


1.92 


1.51 


2.06 


1.08 


2.90 


0.55 


0.98 


0.84 


3.5 





1.75 


1.59 


1.45 


1.01 


0.09 


0.34 


-0.44 


0.92 


-0.67 


4.0 





1.91 


1.70 


1.58 


1.25 


0.18 


0.50 


-0.33 


1.07 


-0.75 


4.5 





1.75 


1.59 


1.59 


1.34 


0.22 


0.58 


-0.25 


1.12 


-0.76 


5.0 





1.78 


1.74 


1.67 


1.35 


0.20 


0.63 


-0.32 


1.15 


-0.72 


5.5 





1.81 


1.95 


1.84 


1.54 


0.29 


0.82 


-0.30 


1.25 


-0.72 




0.0001 


1.80 


1.94 


1.76 


1.60 


0.50 


2.49 


-0.16 


1.10 


0.89 


6.0 





1.76 


1.99 


1.93 


1.72 


0.42 


0.96 


-0.21 


1.30 


-0.76 



Table 5: Elemental abundances [X/Fe] for C, N, O, F, and Ne at the stellar surface at the end of the computed evolution 
for our models with a scaled solar initial composition from Asplund et al. (2009) and Mmix=0.002 Mq for Mass < 3 
and Mmix=0 for Mass > 3 Mr,,. 



Mass (Mq) 


[C/Fe] 


[N/Fe] 


Stromlo 
[0/Fe] 


[F/Fe] 


[Ne/Fe] 


[C/Fe] 


[N/Fe] 


STARS 

[0/Fe] 


[F/Fe] 


[Ne/Fe] 


0.9 


2.74 


1.82 


1.07 


2.78 


2.06 












1.0 


1.86 


0.87 


0.42 


1.52 


0.68 


2.56 


1.37 


1.15 


2.75 


1.88 


1.25 


2.60 


0.98 


0.91 


2.53 


1.84 


2.75 


1.06 


1.22 


2.82 


1.90 


1.5 


2.90 


1.10 


1.17 


3.02 


2.30 


2.84 


0.86 


1.28 


3.04 


2.09 


1.75 


3.10 


1.08 


1.27 


3.40 


2.61 


2.98 


0.83 


1.31 


3.39 


2.35 


1.9 


3.15 


1.14 


1.29 


3.54 


2.72 












2.0" 


3.21 


1.14 


1.31 


3.63 


2.82 


2.78 


0.66 


1.05 


3.27 


2.17 


2.0^^ 


3.05 


1.12 


1.17 


3.38 


2.56 












2.25 


3.27 


1.06 


1.34 


3.73 


2.93 


2.65 


0.51 


0.92 


2.65 


2.05 


2.5 


3.20 


0.85 


1.30 


3.43 


2.77 


2.60 


0.40 


0.84 


3.10 


1.95 


2.5'' 


3.00 


0.80 


1.11 


3.19 


2.49 












3.0 


3.03 


1.46 


0.98 


2.34 


2.00 


2.41 


1.56 


0.65 


1.72 


0.83 


3.5 


2.60 


3.26 


0.95 


1.60 


1.75 












4.0 


2.69 


3.43 


1.29 


2.21 


1.95 












4.5 


2.55 


3.35 


1.20 


1.80 


1.75 












5.0 


2.40 


3.38 


0.91 


0.83 


1.70 












5.5 


2.46 


3.42 


1.00 


0.09 


1.67 












6.0 


2.28 


3.49 


1.24 


-0.45 


1.65 













"For comparison, the M=2 M© model bv ICristallo et all fl2009bf ) produces [C/Fe]=3. 00, [N/Fe]=0.86, [O/Fe]= 1.01, 
[F/Fe]=2.90, and [Ne/Fe]=2.27. ''Model computed with low-temperature opacities from lLederer fc Aringerl (j2009[ ). 



Table 6: Selected key elemental abundances [X/Fe] at the stellar su rface at the en d of th e compute d evolution for 



differen t models. (SM) = Stromlo; (S) = stars models; (BIO) = iBisterzo et all (l2010h : (C09) = ICristallo et al. 



(l2009b[). 



Mass (Mq) 




[Na/Fe] 


[Mg/Fe] 


[Sr/Fe] 


[Ba/Fe] 


[Eu/Fe] 


[Pb/Fe] 


[Ba/Sr] 


[Ba/Eu] 


[Pb/Ba] 


1.25 (SM) 


0.002 


1.68 


0.50 


1.51 


2.17 


1.08 


2.76 


0.66 


1.09 


0.59 


1.25 (S) 


0.002 


1.62 


0.51 


1.65 


1.59 


0.69 


2.75 


-0.06 


0.90 


1.16 


1.3 (BIO) 


-0.001, ST*^ 


0.66 


0.38 


0.32 


0.78 


0.21 


3.22 


0.46 


0.57 


2.44 




-0.001, ST/12 


0.48 


0.31 


0.92 


2.05 


1.17 


2.85 


1.13 


0.88 


0.80 


1.5 (SM) 


0.0006 


1.78 


0.55 


1.59 


2.38 


1.37 


2.88 


0.79 


1.01 


0.50 




0.002 


2.30 


0.91 


1.64 


2.45 


1.47 


3.15 


0.81 


0.98 


0.70 


1.5 (S) 


0.0006 


1.31 


0.42 


1.46 


1.38 


0.47 


2.49 


-0.08 


0.91 


1.11 




0.002 


1.86 


0.76 


1.64 


1.78 


0.85 


2.91 


0.14 


0.93 


1.13 


1.5 (BIO) 


-0.001, ST 


2.26 


1.69 


1.34 


2.04 


1.28 


4.06 


0.70 


0.76 


2.02 




-0.001, ST/12 


2.27 


1.47 


2.26 


2.73 


1.71 


3.11 


0.47 


1.02 


0.38 


2.0 (SM) 


0.0006 


2.81 


1.45 


1.51 


1.89 


0.98 


2.92 


0.38 


0.91 


1.03 




0.002 


3.02 


1.62 


1.78 


2.42 


1.51 


3.24 


0.64 


0.91 


0.82 




0.002^ 


2.46 


1.24 


1.72 


2.46 


1.40 


3.16 


0.74 


1.06 


0.70 


2.0 (S) 


0.0006 


1.90 


1.09 


0.90 


1.52 


0.67 


2.69 


0.62 


0.85 


1.17 




0.002 


2.26 


1.32 


1.34 


2.04 


1.12 


3.03 


0.70 


0.92 


0.99 


2.0 (BIO) 


-0.001, ST 


1.89 


1.44 


1.34 


1.87 


1.05 


4.12 


0.53 


0.82 


2.25 




-0.001, ST/12 


1.87 


1.19 


2.15 


2.96 


1.98 


3.44 


0.81 


0.98 


0.48 


2.0 (C09) 


-o.oor 


1.77 


1.30 


1.01 


1.55 


0.77 


2.88 


0.54 


0.78 


1.33 



"For the models calculated by IBisterzo et al.l (l2010l ) we also need to indicate the choice of the free parameter used in 
those models to describe the efficiency of the ^^C pocket, where ST corresponds to the standard case, as described in 
detail by those authors. This free parameter accounts for a possible range of ^^C and ^^N profiles in the ^^C pocket, 
resulting in a range of total integrated neutron fiuxes. Thus, it is the main free param eter determining t h e fina l 
s-process abundance distribu tion in these raodels. ^Model computed with the opacities by iLederer fc Aringerl (120091 ). 
'^In the models calculated by lCristallo et al.l (l2009bl ) the ^'^C pocket is included self-consistently using an overshooting 
mechanism with the free parameter /3=0.1, as described in details in that paper. This results in a mixing mass of the 
order of — 0.001 Mq, but decreasing with the TP number, while in our models it is kept constant. Furthermore, the 
resulting proton profile is not linear in logarithmic scale (see Fig. 4 of lCristallo et al.ll2009bl ). as imposed in our models. 



